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We study the effects of shear and density inhomogeneities in the formation of naked singularities 
in spherically symmetric dust space-times. We find that in general neither of these physical features 
alone uniquely specifies the end state of the gravitational collapse. We do this by (i) showing 
that, for open sets of initial data, the same initial shear (or initial density contrast) can give rise 
to both naked and covered solutions. In particular this can happen for zero initial shear or zero 
initial density contrast; (ii) demonstrating that both shear and density contrast are invariant under 
a one parameter set of linear transformations acting on the initial data set and (iii) showing that 
asymptotically (near the singularities) one cannot in general establish a direct relationship between 
the rate of change of shear (or density contrast) and the nature of the singularities. However, one 
can uniquely determine the nature of the singularity if both the initial shear and initial density 
contrast are known. 

These results are important in understanding the effects of the initial physical state and in par- 
ticular the role of shear, in determiming the end state of the gravitational collapse. 
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I. INTRODUCTION 

An important outcome of general relativity has been to show that subject to a number of physically reasonable 
assumptions, the final state of gravitational collapse will be singular in a range of settings (see e.g. [1]). An outstanding 
question concerns the nature of the resulting singularities and in particular whether and under what conditions they 
may be naked or black holes [2,3]. Over the recent years a great deal of effort has gone into the study of these 
questions. These have mainly concentrated on collapse in spherically symmetric settings and involve the study of 
collapse of scalar field [4,5] as well as other matter sources including dust [6,7], perfect fluids [8,9], imperfect fluids 
[10-13] and null strange quark fluids [14,15]. 

All these studies show that the end state of the spherical collapse can be either a black hole or a naked singularity 
depending upon the nature of the initial data and the kinematical properties near the singularity. Apart from the 
mathematical questions concerning the likelihood of each end-state, an important question from a physical point of 
view is how the nature of the end-state of collapse in the spherical settings may depend on physical characteristics 
of the regime under consideration. In particular, we ask how does the inhomogeneity and anisotropy influence the 
process and the final outcome of collapse? 

This issue has been studied by a number of authors. Penrose [16] discusses the possibility that anisotropy might 
play a role in the occurence of naked singularities. Herrera et al. [17] considered a spherically symmetric collapsing 
space-time with an anisotropic fluid (with radial and tangential pressures) and showed that high density contrasts 
increase the radial velocity of collapse. They also proved that the local anisotropy of the pressure might have the 
same effect suggesting that this might lead to naked singularities (depending on some inequalities between radial and 
tangential pressures). 

More recently Joshi, Dadhich and Maartens [18] have also considered this question in the special case of marginally 
bounded Lemaitre- Tolman-Bondi (LTB) collapse and found that the shear associated with the initial distribution of 
matter fully determines the final outcome of collapse in these settings. In particular, they have found that the stronger 
the shear near the singularity the more likely it is to have a naked singularity. Soon after, Chan et al. [19] found a 
self-similar solution with heat-flow which was shear-free and, nevertheless, exhibited naked singularities. 

Here, motivated by these previous results, we make a detailed study of the relation between the nature of the 
end-state of LTB collapse with anisotropy and inhomogeneity, by considering more general settings. We study how 
the choices of shear and density contrast functions can influence the nature of the gravitational singularities. In 
particular, we consider the behaviour of the shear function initially, asymptotically (near the singularity) as well as 
at intermediate times and study the effects of each on determining the end state of gravitational collapse. 
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The structure of the paper is as follows: In section 2 we give a brief description of the spherical dust models. In 
Sections 3 and 4 we consider the effects of the shear on the end state of collapse. In section 5 we consider in turn the 
density contrast. Section 6 contains some remarks about space-time extensions after shell-crossing and initial data 
sets. Finally, we conclude in Section 7. 

II. SPHERICALLY SYMMETRIC DUST SPACE-TIMES 

The inhomogeneous spherically symmetric dust space-time can be represented by the Lemaitre-Tolman-Bondi 
(LTB) line element [20-22] 

,i 

ds 2 = -dt 2 + -^—dr 2 + R 2 (d9 2 + sin 2 6dcf> 2 ), (1) 
1 + E 

where r, 9, <\> are comoving coordinates. The dot and prime denote differentiation with respect to t and r respectively 
and R = R(r,t) and E = E(r) are C 2 real functions such that R{r,t) > and E{r) > — 1. The matter-density is 
given by 

m ' 

P ^ = WW (2) 

where m = m(r) is another C 2 real function such that m(r) > 0. The evolution equation for the case of R < 
(corresponding to gravitational collapse) takes the form 



R=-^+E, (3) 
and can be solved for different values of E in the following parametric forms: 
For E < 0: 

Tfl 

R= 2^E) {1 - C ° Sri) 

(v-sm V ) = ^)l { t c -t) (4) 
m 

where < 77 < 2tt and t c — t c (r) is a third C 2 real function that corresponds to the time of arrival of each shell r to 
the central singularity. Note that here and below, t = t c (r) corresponds to 77 = 0. 



For E > 0: 



771 

R = —(cosh 77- 1), 
2Ei 

(sinh rj - 7/) = — - (t c - t) (5) 



where 77 > 0. 
For E = 0: 



3 



R={^) (fc-t)*- (6) 
In what follows we find it also useful to work with the compact form for the solution to equation (3); 

V 777 777 

where G is a positive real function given by 
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arcsm{s/x) y/l-x 
G(x) = , for 1 > x > 

G(x) = | for x = (8) 



-arcsinh(^J~^x) y/1 — x 
G x = J ^ " . f° r > x > -<*>■ 

Using the coordinate freedom to rescale 

R(0,r)=r, (9) 

at an initial time t = 0, equation (7) gives 

t c (r) = ^M, (10) 



m 

where p — —rE/m. So, the initial data set is given by X — {m(r), E(r)}. 

We note that the metric (1) can be matched at a boundary, say r = = const., to the Schwarzschild metric in 
the exterior region. Thus the scenario here is that of a collapsing compact matter region matched to an exterior 
Schwarzschild space-time. 

We shall refer to a singularity as naked if there is a family of future directed non-spacelike geodesies which terminate 
at the singularity in the past. In spherical symmetry, if there are no radial-null geodesies emerging from a central 
singularity then the singularity is necessarily censored [23] . Therefore we shall only be concerned with outgoing radial 
null geodesies which, as can be seen from (1), correspond to the solution of the differential equation 



dt R' 



dr VI + E' 

One can rewrite this equation as (see [6]) 



(11) 



dR 1 / . -dt\ ( E + A/X\ TTfv . 

lR + R lz) = \ 1 -\ i , L )H(X,u), (12) 



where 



du u' V dr I \ V 1 + E 



H(X,u)=( Vu ~P u )X+ (13) 



e u -(r, u -^f3 u )xiG(-PXyj \JP+^ 



and 



R m! u „ E' 



u 



X = — , f] = r — , r\ u = 77—, (3 = r — , /?„ = /?- , 
u m ru E rw 

P=^,P = -r*,e u = e^ 1 ,A=™, (14) 
m m y/uu' u 

_ Jm ... 1 + 8 — n , 3 _ . . 
e=^t>(r) = ^=JL + (r ] --f3)G {p ), 

with the positive real function u = u(r) being monotonically increasing and such that u(0) = 0. Later on we will 
specify u(r) = r 3 . In the cases where E{r) = we will take /?(r) = 0. Following [6], we use the algebraic equation in 
X (obtained from (12)) 



j ^,0,-^-0, (15) 

where the subscript '0' denotes the limit of the associated functions as r — * (respectively u — > 0) along the geodesic, 
in order to demonstrate the existence of radial-null geodesies emanating from the singularity in spherical symmetric 
dust collapse. 

An important outcome of past studies (e.g. [6]) is that the occurrence of black holes or naked singularities as final 
outcomes of collapse depend on the choice of initial data {m(r), E(r)}. This initial data set has to satisfy certain 
regularity conditions namely: 
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1. 

2. 
3. 



there are no trapped surfaces initially, i.e. R(0,r) > m(r), for r > 0. 
there is no shell-crossing during collapse, i.e. R'(t,r) ^ 0, for r > 0. 
the initial matter-density p(0, r) is non-zero at the center r = 0. 



III. DEPENDENCE OF THE NATURE OF SINGULARITY ON SHEAR a(T,R) 

In this section we study the shear a(t,r), in particular its behaviour in the approach to the singularity as well as 
the effects of changing a(t, r) on the nature of the singularity. 

For the LTB metric the shear function (relative to the fluid congruence given by the 4— velocity vector field) is given 

by 



<j(t,r) 



3 I R R' 



(16) 



and its time derivative by 



R 
R 



R 
R 



R/_ 
R' 



R/_ 
R' 



(17) 



We start by recalling the following well-known result (see e.g. [7]): 

Result 1 Consider a dust spherically symmetric collapsing space-time satisfying the regularity conditions. Ifa(t,r) = 
then the final state of collapse is a covered singularity. 

To proceed, we shall use the parametric solutions of the evolution equation to calculate R/R and R! jR! in cases E < 
and E > 0. The case E = has already been considered in [18]. 



In this case we have 



together with 



A. CASE E < 



R = V-E- 



sm?7 
cos 77 — 1 



(18) 



R' 



2^E)) (1 - C ° S??)+ 2(^) Sln ^ 



(19) 



and 



drj . 



R' = (V=E)' 



sm?7 
cos 77 — 1 



where — is given in the Appendix. We can then write 
or 



1 —E 



1 drj 
cos 77 — 1 dr ' 



(20) 



3 \ m J \ (cos 77 — 1 



ML s i n „ _ *L 

2E blU dr 



!) 2 ' (f -^) (cos r?-l) 2 + sin V (cos r? - l)g 



It can be seen that a(t, r) is invariant under the transformation 



(21) 



{m,E) 1 ^ (a 3/2 m,aE), 



(22) 
(23) 
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for an arbirtrary free parameter a > 0. Note that t c (r) does change under this transformation, and so it is a non- 
trivial mapping on Z, i.e. it leads to a different LTB space-time. So, the same shear function a(t, r) can correspond 
to different initial data functions m and E. 

Since we are mainly interested in the behaviour of a(t, r) around the singularity we use an expansion around 77 = 
to obtain 

2\/3(-£) 3 / 2 ( 6 1 -2{AB)' + \AB' - BA' 7 , 2 



r ) ^ 3^T U + ^ + ab\ 120" ) + °« ^ 



A 771 AT, 2 (-^) i 

where A = — — — and B = 

2(-E) m 



B. CASE E > 



In this case, we have 



R=VE-^l-, (25) 
1 — cosh 77 



together with 

drj 



^^(^'(cosh.-D + ^sinh^, (26) 



and 

sinhry r— 1 dt] 



R> = (VeY-^- - VE—^-^, (27) 
1 — cosh 77 f — cosh 77 ar 

where — — is given in the Appendix. We can then write 
or 

a( „ r) = ^ f - Sinh ^ + i^-l 1 ^ (28) 

3 I ™ / V (C0sh 77-1) 2 (f -^)(1-C0sh77) 2 +Sinh7 ? (l-C0sh7 7 )|2 y / ' 

We note that again a{t,r) is invariant under the transformation $0(777, E) = {a 3 ^ 2 m,aE), for an arbitrary free 
parameter a > 0. As a result, the same shear function a(t, r) can correspond to different initial data functions m and 
E. We use the above expressions and their expansions around 77 = to get 

. , 2y/ZE'i ( 6 1 — 2(AB)' - BA' 7r\ \ 2n 

A 771 A D 2El 

where A = —— and B 
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C. Discussion 



The above results show that the shear a(t, r) for both E > and E < cases is invariant under the transformation 
<& a (m,E) = (a 3 / 2 m,aE), with a > 0. This is not true for 6 U . The fact that 9 M is not invariant under $ a has the 
important consequence of allowing the same shear function to produce both naked singularities and black holes. As 
a simple example we take m(r) = m 3 r 3 + mar 6 and E(r) = E 2 r 2 + E 5 r 5 . In this case, letting u — r 3 we obtain from 
(15) the quartic equation 

2A 4 + ^X 3 - 3Q U X + 3V^6„ = 0, (30) 

with 
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e« = 



E5/E2 - m e /m 3 



Pa 



me 
m 3 



3E5 
2E 2 



G( P0 ). 



(31) 



Taking m 3 = 0.1, m 6 = — 0.9, £2 = —0.05, E 5 = 1 we obtain positive roots for (30) which correspond to naked 
singularities. Applying $ a with a = 0.61 to the previous initial data, we get black hole solutions instead and yet the 
same shear function a(t,r). Thus shear a(t, r) does not uniquely determine the final outcome of collapse. 
In the approach to the singularity, the shear in both E > and E < cases can be approximated by 



a(t,r) = - 



V32|£|i 



£)+°fo-'>- 



(32) 



This demonstrates that the behaviour of a(t, r) near the singularity for E ^ depends in general on both m and E, 
whereas in the E = case it depends only on m. 

We can now study the shear strenght in the approach to the singularity depending upon the values of \E\? jm as 
well as the relationship between \E\?/m and the time of horizon formation tn- The surface t — tn(r) is defined by 
R = m. Using (7), this gives 



t H =t c - mG(-E). 



Using the expression (10) for t c , we have 



1 r J 



dt H = 

dm 2 m 3 / 2 



(G(p) + 2pG'(p)) - G(-E), 



where p = —rE/m. Now G is positive, and it turns out that G(p) + 2pG'(p) is also positive, and so 



dm 



< 0. 



(33) 



(34) 



(35) 



Thus for fixed E, regardless of sign, tn does increase with decreasing m. 

It is known that delays in the time of horizon formation can be associated with the increased likelihood of finding 
naked solutions. From (35) we find that smaller m produces delays in tn- Since smaller m also induces stronger 
shear around the singularity (see equation (32)) one may be tempted to relate stronger shears with delays in tn and 
consequently with increased likelihood of finding naked solutions, as was done in [18] for E = 0. However, for E =/= 0, 
this is not necessarily true since in general one needs both m and E to decide whether a solution is naked or not. 

Now recall that for LTB, given an initial function m one can always find a function E such that a naked singularity 
solution exists [6,24]. So we may ask whether smaller m values make it easier to find E functions which give rise to 
naked singularities. In order to answer this question we also note that the part of the singularity surface t c (r) which 
can produce naked singularities is t c (Q), so we are mainly interested in the behaviour of <r(t,r) around r = 0. 

We take C°° functions m and E and integers N > 2, M > 3 so that we can always write m(r) = J2iL 3 m i r % with 
7713 7^ and E(r) = Y^?=2 Ej rJ , with E 2 ^ 0. In this case, we get around r = 



a{t,r) 



V32\E 2 \i 



m 3 



and we find naked singularity solutions from equation (15) if and only if (see e.g. [6,7]) 



Quo € I = 



0, ml 



5V3 N 



13 5^' 



, +00 



(36) 



(37) 



Now a measure for the interval R + \ I is 



5\/3r, 



It is then clear that the smaller m 3 the smaller will be the measure of R+ \ I and consequently, there might exist 
examples where it is easier to find functions E so that Q Uo e I. In fact this is the case for a number of known examples 
(see e.g. [7] for a graphic illustration). However, this phenomena is not general. The problem arises from the fact 
that 9„ depends not only on m 3 and E 2 but also on the higher order coeficients mi and Ej which although do not 
appear in the asymptotic (r — > 0) expression for <r(t,r), can nevertheless be crucial to determine whether Q Uo e I. 
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This then tells us that the information in the vicinity of the singularity given by E2 , m.3 , which determines the shear 
rate in (36), is not enough to determine the nature of the singularity. It is therefore crucial to be sufficiently far from 
the singular region in order to pick up information about higher order " inhomogeneous" coeficients which characterise 
the initial data. 

We note that Joshi et. al. [18] have calculated a(t, r) and cx(0, r) for the case E = and concluded that in approach 
to the singularity these two quantities change with the same rate in r. For the case E 7^ 0, however, one cannot obtain 
expressions for a(t,r) which explicitly depend on r and t so one cannot compare analytically the rates of change of 
cr(t,r) and <x(0,r) in the approach to the singularity (as was done for E = 0). Nevertheless, we prove in the next 
section that the conclusions for E = do not necessarily carry over to the general {E 7^ 0) case. 

We also note that although the function R has the same asymptotic behaviour in the approach to the singularity for 
all E, the shear function a(t, r) docs not, since the dependence of a(t,r) on the function E persists all the way down 
to the singularity for the cases where E 7^ 0. Furthermore, although R "forgets" the E dependence in the approach 
to the singularity, the geodesic equations do not and this is crucial in the study of the nature of the singularities. 

To summarize, in this section we have proved that (i) o~(t,r) is invariant under the linear transformations $ a , (ii) 
a(t, r) diverges in the approach to the singularity at a rate depending on \E\ 3 / 2 /m and (iii) a(t, r) does not uniquely 
specify the nature of the singularity. We have also compared our results for E 7^ to the ones obtained for E = by 
Joshi et al. [18]. 



IV. DEPENDENCE OF THE NATURE OF SINGULARITY ON INITIAL SHEAR a(0,R) 



In this section, we study the initial shear cr(0, r) as well as its possible influence on the nature of the singularities. 
For t = 0, given i?(0, r) = r, one can derive 



<j(0,r) 



3 V r 



E' E 3m w! 



The initial radial velocity of collapse is Vi = R(0, r). So, one can write 

V3 



(38) 



(39) 



which implies that <r(0, r) increases with «j. This confirms, for the case of dust, the results of [17] according to which 
local anisotropy increases the radial velocity of collapse. Now, for R(t, r) 7^ 0, we have that tr(0, r) = if and only if 



2 r 2r 2 2 r 



(40) 



Result 2 Consider a dust spherically symmetric space-time (M, g) with ct(0, r) = 0. Then the following are equivalent: 

1. <r(0,r) = 0. 

2. (M, g) is spatially homogeneous and isotropic. 

3. a(t,r) = 0, for all t > 0. 

Proof: Using (17) for <r(0, r) = we obtain 

(41) 



(42) 



implies E(r) = E 2 r 2 . These correspond to spatially homogeneous initial data (see e.g. [7]). On the other hand, spatial 
homogeneity implies &(t, r) = • 





&(t,r) = 




(i 


-i) 










which gives 












a(0,r) - 


6 


ml 


3m 




Note that this does not depend on E. 


The condition 


<t(0, 




implies m(r) = m^r 3 
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Result 2 indicates that, as expected in general, even if we start with zero initial shear in an inhomogencous space-time, 
shear will be generated after evolution starts. This can also be seen from the relationship between a(0, r) and p(0, r): 

(43) 

This, however, is not the case for E = 0: 

Result 3 Consider a dust spherically symmetric space-time with E = 0. Ifa(Q,r) = then a(t,r) = 0. 

Proof: For E = 0, if <r(0, r) = then one must have m = m 3 r 3 which gives a(0, r) = 0. Now apply Result (2) • 

For the case of E = 0, since shear depends only on m there is a unique correspondence between m and tr(0, r) 
and one can therefore substitute the initial free function m by er(0, r). One can then study how the initial shear influ- 
ences the final outcome of collapse, as in [18]. In particular, if er(0, r) = the final outcome of collapse is necessarily 
a black hole. An interesting question is whether this result necessarily holds for £ ^ 0. 

The situation is in fact more complicated in the case E ^ 0, where as we shall show the same initial shear can give 
rise to different outcomes. Before doing so we note that one can fine tune the initial data so that c(0,r) = and 
&(0,r)^0. 

Result 4 Consider a dust spherically symmetric collapsing space-time with m(r) = r 3 g(r) and E(r) = r 2 Q(r), such 
that g and Q are C 2 real functions in [0, rj] with g(0) ^ and Q(0) ^ 0. Then a(0, r) = if and only if Q 1 + g' = 0. 

Proof: Follows directly from expression (38) • 



<r(0,r) = 



V3 



p(0,r)- 



3m 



Therefore there can exist situations where the initial shear is zero and nevertheless a naked singularity will be formed. 
Consider, for example, an initial data set T given by the functions m(r) = m 3 r 3 + m§r 6 and E(r) = E 2 r 2 + E$r 5 . 
Taking the initial coefficients to be E 2 — m 3 = 0.1 and E 5 = — m 6 = 1 we ensure that c(0,r) = and <r(0,r) 7^ 
(this initial data also satisfies the regularity conditions). Solving (30) for these values one finds a positive root Xq 
which therefore corresponds to a naked singularity. 

Taking the same initial data except for the value of m 3 — 0.2 we find a black hole, while we still have cr(0,r) = 
and &(0, r) 7^ 0. One can also find examples for E 2 < such as E 2 = — 0.05,m 3 = 0.1, £5 = -m 6 = 1 giving a 
naked singularity and E 2 = —0.05, m 3 = 1, E5 = — m 6 = 1 giving a black hole, both initial sets with ct(0, r) = and 
*(0,r)^0. 

Furthermore one can find open intervals in the coefficients of m and E such that the initial shear is zero and 
yet naked singularities form as in the next example: Taking m(r) = m^r 3 + mar 6 and E(r) = E 2 r 2 + E 5 r 5 with 
me + E 5 = 0, then from (31) we obtain 

e„ = J (V^o + iPo |)G(po)) , (44) 

with po = —E 2 /m^,. Taking the elliptic case for which < po < 1, i.e. 

£ 2 e[-m 3 ,0[, (45) 

and consequently 2/3 < G(po) < n/2. In this case y/1 — po + (po — 3/2)G(po) is always bounded and negative. So, in 
order to have Q uo > and therefore avoid shell-crossing we take E 5 > 0. Furthermore in order to ensure that Q ug 

falls in one of the intervals given by (37), e.g. O M0 > m^fy + ^p), we impose 



E 5 G 



E 2 m 2 3 /3 (f + Sffi 



yr^+(po-3/2)G(p ) 



. +00 



(46) 



So, cr(0, r) = and naked singularities as final state of collapse can be achieved for open subsets (in the sense made 
precise above) of the initial data. The possibility of having the same initial shear with different end states for the 
collapse can also be found for non-zero initial shear. This can be shown by generalising Result 4: 

Result 5 Consider a dust spherically symmetric collapsing space-time with initial data as given in Result 4- Then 
two different sets of initial data {m^, E^} and {m/ 2 -*, E^} have the same initial shear function <i(0,r) if and only 
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An example of this is the naked singularity solution obtained for i = 6 and j = 5 with 777,3 = 0.1, E 2 = —0.05, mg = 
— 0.9, E 5 = 1, and the black hole solution 7773 = 1,E 2 = — 0.95,m 6 = — 0.9, £5 = 1 both having the same non-zero 
initial shear profile. Again one can find open subsets of the initial data which have the same non-zero initial shear 
and yet different outcomes. 

To summarize, in this section we have proved that (i) the same initial shear can result in two different outcomes 
black holes and naked singularities depending on {m, E}, so initial shear does not uniquely determine the final state 
of collapse; (ii) it is possible to have zero initial shear and a naked singularity endstate, and (iii) these results can be 
obtained for open sets of initial data. 



V. DEPENDENCE OF THE NATURE OF SINGULARITY ON INHOMOGENEITY 



The studies of the previous sections for the shear can now be repeated for inhomogeneity measures. In this section, 
we study the influence of the density contrast on the nature of the singular endstates. We shall use the pointwise 
density contrast measure 

777" rn'i?" 2m' 
9 ( *' r) _ 1VR 2 ~ R 2 (R') 2 ~ W ( ' 

and start by recalling the following well-known result: 

Result 6 Consider a dust spherically symmetric collapsing space-time satisfying the regularity conditions. Ifp'(t,r) = 
then the final state of collapse is a black hole. 

We shall therefore consider the cases where p'(t, r) 7^ 0. If we now use the formulae in the Appendix and substitute 
for R,R' and R" in equation (47) we obtain, after expanding around 77 = 0, 

P %r) = -m) 3 ^ + 0(^). (48) 
\ 777 J rf 

Thus as in the case of shear, the inhomogeneity measure p'(t, r) depends asymptotically on both E and 777. In order 
to prove that m' (A\E\/m) 3 does not uniquely determine the nature of the singularity we note that this quantity is 
invariant under the linear transformation $ (m, E) = (a 3 / 2 m,aE), for an arbitary free parameter a > 0. In fact, 
by substituting R, R' and R" in (47) we also find that the all time function p'(t,r) is itself invariant under $ . 
Interestingly $ also leaves invariant the shear function <r(t, r) and the examples of section III C can also be used here 
in order to show that the same function p'(t, r) can give rise to both naked singularities and black holes. 
We now study the initial density contrast function p'(0, r). From (47) we obtain for t — and R(0, r) = r 

//„ \ 777" 2m' , . 

?(M = — -— ( 49 ) 

so p'(0,r) does not depend on E, which makes it clear that a given p'(0,r) can correspond to both black holes and 
naked singularities, depending on the choice of E (see e.g [7]). So p'(0, r) cannot uniquely determine the final outcome 
of collapse, except in the E — case. An interesting question is whether there are naked singularity solutions for 
p'(0, r) = 0. In order to answer this question we start by recalling a simple result: 

Result 7 Consider a dust spherically symmetric space-time with p'(0,r) = 0. Then the space-time is spatially 
homogeneous and isotropic if and only if <r(0,r) = 0. 

Proof: If p'(0, r) = then 777 = m 3 r 3 , 7773 € R \ 0, in which case cr(0, r) = if and only if E = E 2 r 2 • 

We recall that in the E = case the space-time is spatially homogeneous and isotropic if and only if p'(0, r) = 0. We 
shall then take the inhomogeneous E ^ cases where it might be possible, for a(0,r) ^ 0, to find initial data such 
that p'(0, r) — and nevertheless have a naked solution. As an example of this we take m(r) = m 3 r 3 , 7773 e R + , and 
E(r) = E 2 r 2 + E 5 r 5 such that 

e -- = f b^-iM ' (50) 

By taking the elliptic case with E 2 € R~ and 
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we ensure both that 9 Uo > and that equation (15) has positive roots corresponding to naked singularities. 

To summarize, in this section we have shown that (i) p'(t,r) is invariant under the linear transformations $ a , (ii) 
p'(t, r) diverges in the approach to the singularity at a rate depending on m'(\E\/m) 3 , (iii) p'(t, r) docs not uniquely 
specify the nature of the singularity, (iv) the same initial inhomogeneity function p'(0,r) can result in two different 
outcomes black holes and naked singularities depending on {m,E}, so p'(0,r) does not uniquely determine the final 
state of collapse and (v) it is possible to have zero initial density contrast and a naked singularity endstate. These 
results can be obtained for open sets of initial data. 



VI. NON-UNIQUENESS, SHELL-CROSSING AND THE INITIAL DATA SET 

Finally, we consider two important issues in this section: (I) Non-uniqueness of space-time extensions after shell- 
crossing singularities and (II) Initial data sets based on initial shear and initial density contrast functions. 

(I) The theme of our discussion has been how certain quantities which have been suspected of determining the 
outcome (naked singularity or black hole) of spherical dust collapse in fact fail to do so. We have emphasised the 
significance of the well-known fact that one needs to know both E and to in order to make this prediction. However 
there is another situation in which even the knowledge of both these functions is not enough to completely determine 
the space-time; that is when the initial data are such that shell-crossing singularities can occur. It has been shown 
recently how one can construct a dynamical extension - i.e. one based on the field equations - through a shell-crossing 
singularity and obtain global weak solutions of the field equations for R > 0, t > [25]. However these weak solutions 
are not unique to the future of the shell-crossing. On the other hand, the earliest point of the shell-crossing is always 
globally naked, regardless of how the extension is constructed, so the non-uniqueness is not so severe as to allow the 
same initial data give rise to either a black hole or a naked singularity. 

(II) As we have seen, neither <r(0, r) nor p'(0, r) alone can uniquely determine the final outcome of collapse. How- 
ever, both functions might form a perfectly good initial data set. The question is how from a set of observables 
J = {p'(0, r), <t(0, r)}, given on the onset of collapse, can one predict (uniquely) the final outcome of collapse. 

Now m(r) is uniquely determined by p'(0,r) through the linear ODE (49): 

to" - ^m' = r 2 p'{0,r), (52) 

with initial conditions m(0) = m'(0) = 0. The shear <r(0, r) can be written as a function of a new variable w(r) = 
E + m/r: 

x v3 i . w 1 ,. . . 

cT(0,r) = ^-w-*(-- + -w'), (53) 

with initial conditions w(0) = w'(0) = 0. So, from <r(0,r) one can determine w(r) by solving the linear ODE 

v'-\v = V3a(0,r), (54) 

where w = y 2 . So, using (52), we can get m(r) from p'(0,r). Knowing m(r) and a(0,r) we can determine E from 
(54). So, if we know both the density contrast p'(0, r) and the anisotropy measure cr(0, r) we can recover both E and 
to, and so predict the outcome of collapse. Since this process relies on solving linear ODEs, then one can ensure that 
solutions exist and are unique, for the given set of initial conditions. 

It is worthwhile noting that even though the set {u(0, r), p'(0, r)} can specify the initial data uniquely, the set 
{a(t, r), p'(t, r)} specified at any other time cannot. The reason is that the latter is invariant under the one-parameter 
set of linear transformations $ a acting on the initial data set {to, E}, while the former is not and this is because even 
though initially one can use the rescaling to define R(0, r) = r, this cannot be repeated at subsequent times. 
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VII. CONCLUSIONS 



We have considered LTB collapsing models in order to study the influence of inhomogeneity and anisotropy in the 
nature of the resulting singularities. We have studied the shear and density contrast in three ways: (i) by considering 
the all time functions cr(t,r) and p'(t,r), (ii) by considering the asymptotic functions a(t, r) and p'(t,r) and (iii) by 
considering the initial functions er(0,r) and p'(0,r). 

We have found that both shear <x(i, r) and density contrast p'(t, r) are invariant under a one parameter set of linear 
transformations acting on the initial data set {to, E}. We have also found that asymptotically one cannot establish, in 
general, a direct link between the strenght of shear (or density constrast) and the nature of the singularities. Finally 
we have found that the same initial shear (or initial density contrast) can give rise to both naked singularities and 
black holes, depending on the choice of the initial density contrast (or initial shear, respectively). In particular this 
can happen for zero initial shear or zero initial density contrast. These results were obtained for open sets of initial 
data. We conclude that neither anisotropy nor inhomogeneity features as given by the shear or density contrast can 
alone uniquely specify the end state of collapse for E ^ 0, in contrast to the results established for E = in [18]. 

However, we have proved that if we know both the initial shear er(0,r) and the initial density contrast p'(0,r) 
functions then, given the appropriate initial conditions, we can determine uniquely the initial LTB data functions 
{to, E} and hence determine the nature of singularity. In this sense one needs information about both the inhomo- 
geneity and anisotropy at the initial hypersurface in order to predict the outcome of the collapse. This is important 
for understanding the effects of the initial physical state in determiming the end state of gravitational collapse. 
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